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Abstract. Let M be a compact connected oriented n— 1 dimensional manifold 
without boundary. In this work, shape space is the orbifold of unparametrized 
immersions from M to R". The results of [T], where mean curvature weighted 
metrics were studied, suggest incorporating GauB curvature weights in the 
definition of the metric. This leads us to study metrics on shape space that 
are induced by metrics on the space of immersions of the form 



Here / e Imm(M,K") is an immersion of M into R" and h,k £ C°°{M,R") 
are tangent vectors at /. g is the standard metric on R", f*g is the induced 
metric on M, vol(/*g) is the induced volume density and $ is a suitable 
smooth function depending on the mean curvature and Gaufi curvature. For 
these metrics we compute the geodesic equations both on the space of immer- 
sions and on shape space and the conserved momenta arising from the obvious 
symmetries. Numerical experiments illustrate the behavior of these metrics. 



Riemannian metrics on shape space are of great interest for many procedures 
in data analysis and computer vision. They lead to geodesies, to curvature and 
diffusion. Eventually one also needs statistics on shape space. The shape space 
used in this paper is the orbifold i?i(Af, E") of unparametrized immersions of a 
compact oriented n — 1-dimensional manifold M into K", see [S]: 



The shape space of submanifolds in M" of type M used in [1] [5] is a special case 
of this. 

The quest for suitable Riemannian metrics was started by Michor and Mumford 
in O E], where it was found that the simplest and most natural such metric 
induces vanishing geodesic distance on shape space Bi. One attempt to strengthen 
the metric is by adding weights depending on the curvature. For curves, where there 
is only one notion of curvature, this was done in |12l I13| . In higher dimensions, 
there are several curvature invariants. We are especially interested in the case 
n = 3, where the curvature invariants are the mean and GauB curvature. Mean 
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1. Introduction 



B,(M,M") = Imm(M,IR")/Difr(M). 
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curvature weights were used in [1 . In this paper we take up these investigations 
by adding a Gaufi curvature term to the definition of the metric. 

A Riemannian metric on Imm(Af, M") is a family of positive definite inner prod- 
ucts Gf{h,k) where / e Imm(M,M") and h,k e C°°(M,]R") ^ r/Imm(M,M") 
represent vector fields on M" along /. We require that our metrics will be invariant 
under the action of Diff(M), hence the quotient map resulting from dividing by 
this action will be a Riemannian submersion. This means that the tangent map of 
the quotient map Imm(M, M") — > i?i(M, E") is a metric quotient mapping between 
all tangent spaces. Thus we will get Riemannian metrics on Bi. Riemannian sub- 
mersions have a very nice effect on geodesies: the geodesies on the quotient space 
Bi are exactly the images of the horizontal geodesies on the top space Imm. 

The metrics we will look at are of the form: 

Gf{h,k)^ f ^{TT{L),det{L))g{h,k)vo\{rg) 
J hi 

Here g is the standard metric on M", f*g is the induced metric on A/, vol(/*g) is 
the induced volume density, L is the Weingarten mapping of / and $ : — )• ]R>o 
is a suitable smooth function depending on the mean curvature Tr(L) and Gaufi 
curvature det(I/). Since <i> is just multiplication with a smooth positive function, 
the horizontal bundle consists exactly of those vector fields which are pointwise 
normal to f{M), see p^. 

Big parts of this paper can be found in Martin Bauer's Ph.D. thesis [3]. 

2. Differential geometry of surfaces and notation 

We use the notation of [TU]. Some of the definitions can also be found in [7]. 
Similar expositions in the same notation are [TJ |^ . 

In all of this chapter, let / : M — )■ Imm(M, M") be a smooth path of immersions. 
By convenient calculus / can equivalently be seen as a smooth mapping / : 
R X M ^ M" such that f{t, ■) is an immersion for each t. We will deal with the 
bundles 

T;M T^M (g) f{t,-)*TR" pTt,T^M (g) f*TW 

M M Rx M 

Here TJM denotes the bundle of ( I )-tensors on AI, i.e. 

r s 

TIM = TM (S) T*M and pr2 : M x Af ^ M. 

Let g denote the Euclidean metric on M". The metric induced on M by f{t, •) is 
the puUback metric g = f{t, )*g. On tensor spaces TJAf we consider the product 
metric gl := g ® g^^ . A useful fact is (see [2]) 



(1) gl{B,C) :^TT:{g-^Bg-'^C) for C £ r°M if B or C is symmetric. 
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Furthermore g^^g yields a metric on TJM (g) f{t,-)*TR'\ Let V denote the 
Levi-Civita covariant derivative acting on sections of all of the above bundles. The 
covariant derivative admits an adjoint, which is given by V* = — Tr^ V. For more 
details see 

The normal bundle of an immersion f(t, •) or a path of immersions /, respectively, 
is given by 

or 

j Y 

M Rx M 

Let v designate the positively oriented unit normal field with respect to the 
orientations on M and M". Any field h along f{t,-) or / can be decomposed 
uniquely as 

h = Tf{t,-).h^ + h-^M 

into parts which are tangential and normal to /. The two parts are defined by the 
relations 

= g{h,iy) e C°^{M) 

e X(Af), such that g{h^ ,X) = g{h,Tf{t,-).X) for all X e X{M). 

Let denote the second fundamental form of f{t, •) on M and L the Weingarten 
mapping, which are related to the Levi-Civita covariant derivative V on (M",^) by 
the following equation: 

g(Vxiy,Tf{t,-).Y) = ~siX,Y) = ~g{LX,Y). 

The eigenvalues of L are called principal curvatures and the eigenvectors principal 
curvature directions. Tr(L) is the mean curvature, and in dimension two det(L) is 
the Gaufi-curvature (we use this name generally). 

3. VARIATIONAL FORMULAS 

We will calculate the derivative of the GauJ3-curvature. The differential calculus 
used is convenient calculus as in [5] . Proofs of some of the formulas in this chapter 
can be found in [1 HH H H] and in [H] . 

3.1. Setting for first variations. Let us consider a function F defined on the 
set of immersions Imm(M, M"), an immersion /o and a tangent vector h to /q. 
Furthermore we choose a curve of immersions 

/ : M ^ Imni(A-f,E") with /(O) = fo and dt\of = h = h^.v^" + Tfo.h^. 

In this setting, the first variation of F is 

DuoM)F^Tf„F.h = dt\oF{f{t)). 
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3.2. Tangential variation of equivariant tensor fields. Let the smooth map- 
ping F : Imm(A/, M") r(TJM) take values in some space of tensor fields over 
M, or more generally in any natural bundle over M, see [S]. If _F is equivariant 
with respect to puUbacks by difFeomorphisms of M, i.e. 

F{f)^{rF){f)^cf>*{F{{cb-'rf)) 

for all (j) G Diff(Af) and / € Imm(M, M"), then the tangential variation of F is 
^(/o,T/o.hT)F = dtloFifo o Flf) = dt\oF{{Flfrfo) 
^dMFlfnF{fo))^C,TiF{fo)). 

3.3. Variation of the volumeform [1] section 3.5]. The differential of the volume 
form 

( Imm ^ r2"-i(Af), 

[ f ^ vol(ff) = vol(rg) 

is given by 

i?(;„,;,)Vol(/*g) = 9t|ovol(r5) = (div3°(/iT)-Tr(L/«)./i^)vol(5o). 

3.4. Variation of the Weingarten map [1] section 3.8]. The differential of the 
Weingarten map 

( Imm ^ r(End(TM)), 
[ f ^ Lf 

is given by 

D^j^^h)Lf = dtloL = go'^'ih^) + h^{Lf"f + C^,r {Lfo). 

3.5. Variation of the Gauficurvature. The differential of the Gaufl curvature 

Imm ^ C°°(M), 
/ ^ det(L^) 

is given by 

D^f^^h) det(L^) = Tr(L^").det(L-^«)./i^ +g^(g.C(L-^°),V2(/i-L)) + dAel{Lf''){h^) 

where C{L) is the classical adjoint of L uniquely determined by 

G{L).L = L.C{L) = det{L).I 

Proof. For the normal part we have 

atlodct(i) =Tr(C(L-^°).at|oL) 

- Tr (c(L-^«).((/*.g)-iv2(;,^) + h^{Lf« f)"^ 

= Tr ( C(L^«).(/o5)"^V^(/i-^) + det(i-^«).L^" " 
= Tr{L^»).Act{LS»).h^ + gl(g.C{Lh),V\h^)) 
For the tangential part, we use the observation of section [3^ □ 



CURVATURE WEIGHTED METRICS ON SHAPE SPACE OF HYPERSURFACES IN n-SPACE5 



4. The geodesic equation on Imni(M,E") 

4.1. The geodesic equation on Inim(M, M"). We derive the geodesic equation 
only for G'*-metrics witli $ = $(det(i)). This equation can be combined with 
the geodesic equation for metrics weighted by mean-curvature and volume [H sec- 
tion 5.1]. The resulting geodesic equation for <I> = $(Vol, Tr(i), det(L)) on Imm is 
very long, but on Bi it is much shorter and we print it in section [5. 2| 



We use the method of [II, section 4] to calculate the geodesic equation as 
ftt = -j^Hf{ft,ft) ~ Kf{ft, ft), 
where H, K are the metric gradients given by 

Df,^Gf{h, k) = GfiKfim, h), k) = G^(m, Hf{k, h)). 

So we need to compute the metric gradients. The calculation also shows the exis- 
tence of the gradients. Let m G Tf Imm(M, K") with m = .v^ + Tf.mJ. 

D^f^,^^Gj{h, k) = / $(det(i)) -g{h, k) vol(<?) 

JM 

$'(det(L))(i^(/,„) det(i)) g{h, k) vol(g) 

M 

+ [ $(det(i))g(/i,/c)p(;,„)V0l(5)). 

JM 

To read off the i^-gradient of the metric, we write this expression as 
^-9 dot(X)) + k) vol(,) 

Therefore, using the formulas from section [3] we can calculate the K gradient: 
XKr.,/.)^[|(i.(,,„)det(X)) + ^i/^^ 

^ ( Tr(L). det(L).m^ + gl{g. C(L), V'^{m^)) + dAei{L){mJ) 

-f divf(rn^) -Tr(L).TO-^ 

To calculate the 77-gradient, we treat the two summands of D(^f „i-jGj {h, k) sepa- 
rately. The first summand is 



M 



det(L))5(/i,fc)vol(5) 



= ( <^'Tv{L).dei{L).m^.g{h,k)Yo\[g)+ [ <i>'ddet{L){m^)g{h,k)vol{g) 

JM JM 

+ [ g°{<P'.g{h,k).g.C{L),\/\m^))vo\{g) 

JM 

= [ <i>'TT{L).det{L).m^.g{h,k)vol{g) 

JM 



M 



$'.g(gradS(det(L)),m^)5(/i,fc)vol(5) 
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+ ^ V* V* ($'.5. C{L)g{h, fc)) .m^ vo\{g) 



= Gj l^m, — Tr(L). det(L).5(/i, k)M + -^g{h, k).Tf. gradS(det(L)) 

+ g.C{L)g{h,k))M 

In the calculation of the second term of the Hf(m, h) gradient, we will make use 
of the following formula from [T] section 5.1], which is valid for € C°°(Af) and 

X e je(M): 

/ dc^[X)Yo\{g)^- ( 0.div(X)vol(5). 

JM J hi 

Therefore we can calculate the second summand, which is given by 

$.g(/i,fc)py,^)V0l(5)) = / $.g(/i,A:)(divS(m^)-Tr(i).m^)vol(g) 

M J M 

d{<i>g{h, /fc))(m^) vol(g) + Gj{m, -g{h, k) Tr(i).i/) 



M 



M 



g{Tf. gY&d^<f g{h, k)),m) vol(g) + Gf{m, -g{h, k) Tr(X).z.) 



G* (™, grad«(<i>ff(/i, fc)) - fc) Tr(L). 



Summing up all the terms the if-gradient is given by 



Hf{h,k) 



^ . Tr(L). det(L).5(/i, fc) + - V* V* ($'.5. C{L)-g{h, k)) - g{h, k) Tr(L) 



Tf. 



'P'g{h, k) grad^(det(L)) - grad^($.g(/i, fc)) 



The geodesic equation for a Gaufi- curvature weighted metric G* with $ = $(det(L)) 
on Imm(Af, M") is then given by 



ht^\[^- Tr(i). det(L). \\h\\^ + ^ V* V* ($'.5. C(L) H/if) - 
^ ~" V||/i||^grad9(det(L)) -grad5($. 



Tr(L) 



2$ 



T/. 



[- 



Tr(L).det(L)./i-L + ^^(.g. C(L), V2(/i-L)) +ddet(i)(;i^ 



+ div5(/i^) - Tr{L).h^]h. 



4.2. Momentum mappings. The metric G is invariant under the action of the 
reparametrization group Diff(M) and under the Euclidean motion group K" x 
SO(n). According to [TJ section 4] the momentum mappings for these group actions 
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are constant along any geodesic in Imm(A/, M"): 

vol((7) € r(r*M (g)Af vol(M)) reparam. momentum 
<i>/f vol(.g) e (M")* linear momentum 

^if ^ ft) vo\{g) e /\^M" ^ so(n)* angular momentum 



5. The geodesic equation on i3i(M, M") 

5.1. The horizontal bundle and the metric on the quotient space. Since 
vo\{f *g), Tr{L) and det(_L''') react equivariantly to the action of the group Difl'(il/), 
every G^-metric is DifF(M)-invariant. Thus it induces a Riemannian metric on Bi 
such that the projection tt : Imm Bi is a Riemannian submersion. For every 
almost local metric G* the horizontal bundle at the point / equals the set of sections 
of the normal bundle along /. Therefore the metric on the horizontal bundle is given 

by 

Gj{a.v,b.u)^ ( $.a.6vol(/*5). 
See [U section 6.1] for more details. 





5.2. The geodesic equation on Bi{M, M"). The calculation of the geodesic equa- 
tion can be done on the horizontal bundle instead of on Bi. This is possible because 
for every almost local metric a path in Bi corresponds to exactly one horizontal 
path in Imm (see [T] ) . Therefore geodesies in Bi correspond to horizontal geodesies 
in Imm. A horizontal geodesic / in Imm has ft = a.v^ with a G C°°{R x M). The 
geodesic equation is given by 

ftt = = ^H{a.iy, a.v) - K{a.v, a.v), 

normal tang. 

see [T] section 4]. This equation splits into a normal and a tangential part. The 
normal part is given by 

at = g{-H{a.i', a.v) — K{a.v, a.v), v). 

From the conservation of the reparametrization momentum, see IJ it follows that 
the tangential part of the geodesic equation is satisfied automatically. 

Therefore the geodesic equation on Bi for $ = $(det(L)) is given by 



ft = a.v^, 



1 



I TT{L).a' + iv* V* C{L)a') - <l>'g° (<?• C{L),V^{a)) .a 

- y.Tr(L).det(L).a2 
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Again the geodesic equation for a more general almost local metric G where 
$ = <i>(Vol, Tr(_L), det(L)) can be obtained by combining the above equation with 
the results in [IJ section 5.1]. It reads as 



ft = a.i'^, 






1 


|Tr(L).a2-iTr(L) ljd,<i>)a' 


vo\{g) + {di<^>)a. [ Tr(L).a vol(g) 

JM 




- -A((a2$)a') + (92$) 


.aAa 


- —a Tr(i ) 




+ ^V*V*((a3<i>)..g.C(i) 




{d^^)gl(g.C{L),V\a)).a 




- ^.Tr(L).det(L).a2 







For the case of curves immersed in M^, this formula specializes to the formula 
given in [T3] section 3.4]. (When verifying this, remember that A = — and 
Tr(L) — det(L) — k in the notation of |13j.) 



6. Numerical results 

We want to solve the boundary value problem for geodesies in shape space of 
surfaces in M'^ with respect to curvature weighted metrics, more specifically with 
respect to G^-metrics with $ = 1 + A. Tr(L)2'= + B. det{L)'^K In [T] we did this for 
metrics depending on the mean curvature only. 

We will minimize horizontal path energy 

^'°'(/)= / $(Tr(i),det(L))(/,^)2vol(5) 

Jo JM 

over the set of paths / of immersions with fixed endpoints. By definition, the hor- 
izontal path energy does not depend on reparametrizations of the surface. There- 
fore we can add a penalty term to the above energy that controls regularity of the 
parametrization. 

We reduce this infinite dimensional problem to a finite-dimensional one by ap- 
proximating immersed surfaces by triangular meshes. We solved the resulting op- 
timization problem using the nonlinear solver IPOPT (Interior Point OPTimizer 
[15]). IPOPT was invoked by AMPL (A Modeling Language for Mathematical 
Programming [HI ). The data file containing the definition of the combinatorics 
of the triangle mesh was automatically generated by the computer algebra system 
Mathematica. 



6.1. Discrete path energy. As in T, section 10.1] we define the discrete mean 
curvature at vertex p as 

, , 1 1 vector mean curvature 1 1 || V (surface area) || 
Tr(I/) = — — 



|vector area|| ||V(enclosed volume)|| 
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Here V stands for a discrete gradient, and 

(vector mean curvature)^ = (V(surface area))^ — (cota^ + cot /3i){p — Pi) 

is the vector mean curvature defined by the cotangent formula. In this formula, 
ai and Pi are the angles opposite the edge {p,Pi) in the two adjacent triangles. 
Furthermore, 

(vector area)p — (V(enclosed volume))^ — i^(/).( surface area of /) 

/eStar(p) 

is the vector area at vertex p. The discrete Gauss curvature can be defined as 

'^^^^^^^^^ = Areaoflr(p) - 
Here Q{p) stands for the angular deflection at p, defined by 

#(star(p)) 

where 9i denotes the internal angle of the i-th corner of vertex p and #(star(p) the 
number of faces adjacent to vertex p. Using these definitions, the horizontal path 
energy and the penalty term were discretized as in [U section 10.1]. 



6.2. Geodesies of concentric spheres. In this chapter we will study geodesies 
between concentric spheres for Gaufi curvature weighted metrics, i.e. metrics with 
$ = 1 + i?det(L)^'. For mean curvature weighted metrics this has been done in [T]. 

The set of spheres with common center x E M" is a totally geodesic subspace of 
Bi with the G^-metric (see jlj section 10.3]). Within a set of concentric spheres, 
any sphere is uniquely described by its radius r. Thus the geodesic equation reduces 
to an ordinary differential equation for the radius. It can be read off the geodesic 
equation in section [5] when it is taken into account that all functions are constant 
on each sphere, and that 

rt = a, L^-^ldTM, Tr(L) = -^^, det{L) = 

Then the geodesic equation for $ = 1 + B dct(L)^' on a set of concentric spheres 
in Bi reads as 

_ ^n-l fl 21. B 

Note that det(L) = (Tr(L)/(n — 1))" ^. Therefore this equation is equal to the 
equation for metrics weighted by mean curvature with suitably adapted coefficients 
(see dJ section 10.3]). 

This equation is in accordance with the numerical results obtained by minimizing 
the discrete path energy. As will be seen, the numerics show that the shortest path 
connecting two concentric spheres in fact consists of spheres with the same center, 
and that the above differential equation is satisfied. Furthermore, in our experiments 
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the optimal paths obtained were independent of the initial path used as a starting 
value for the optimization. 

A comparison of the numerical results with the exact analytic solutions can be 
seen in figure[l] The solid lines are the exact solutions. For the numerical solutions, 
50 timesteps and a triangulation with 320 triangles were used. 




Figure 1 . Geodesies between concentric spheres for Gauss curva- 
ture weighted metrics with $ = 1 + Bdet{L)^. 



The space of concentric spheres is geodesically complete with respect to the G 
metrics, $ = 1 + i3det(L)^', if I > 1- To see this, we calculate the length of a 
path / shrinking the unit sphere to zero [n — i): 

Lf (/) = I' ^lGf\ffjf)dt = £ ^ljl + B.detiL)^^)rfvol{g)dt 

= \rt\\l{l + B. det{LY^)A:rHdt = ^Vt^ ry(l + B.-^)dr. 
The last integral diverges if and only if Z > 1 . 



6.3. Translation of a sphere. For metrics with <I> = 1 + B. det(L)^', the following 
behaviours can be observed: 

• For spheres of a certain optimal radius, pure translation is a geodesic. 

• When the initial and final shape is a sphere with a non-optimal radius, the 
geodesic scales the sphere towards the optimal radius. 
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• When too much scaling has to be done in too little time (always relative to 
the parameter B ), then the geodesic passes through an ellipsoid-like shape, 
where the principal axis in the direction of the translation is shorter, see 
figure 







c 






m 




c 











Figure 2. Geodesies of translation for $ = 1 + Bdet{L)'^, with 
B = 1 (first row) and B = 20 (second row) progresses from left to 
right. 



To determine under what conditions pure translation of a sphere is a geodesic, 
let /t = /o + b{t) ■ ei, where /o is a sphere of radius r and where b{t) is constant on 
M. Plugging this into the geodesic equation 4.1 yields an ODE for b{t) and a part 
which has to vanish identically. The latter is given by: 



2 1 



$- = 
r 



For $ = 1 + B. det(L)2', Z > 1 this yields the solution r = %/B{2l ~ 1). 



6.4. Deformation of a sphere. For $ = 1 + i?det(L) small deformations look 
much the same as for conformal or mean curvature weighted metrics ( compare 
[TJ figure 14]/ Larger deformation are somewhat irregular, but adding a mean 
curvature term yields smooth geodesies again (see figure^. 



V C C L' 



Figure 3. Geodesic between two deformed shapes for <i> = 1 
Tr(L)^ + det(_L)^. Time progresses from left to right. 
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